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Abstract 

A major controversy has arisen in QCD as to how to split the total angular momentum into 



i-C ' separate quark and gluon contributions, and as to whether the gluon angular momentum can itself 

be split, in a gauge-invariant way, into a spin and orbital part. Several authors have proposed 

<u : . 

various answers to these questions and offered a variety of different expressions for the relevant 



operators. I argue that none of these is acceptable and suggest that the canonical expression for 
the momentum and angular momentum operators is the correct and physically meaningful one. It 



> 

l/V 

, is then an inescapable fact that the gluon angular momentum operator cannot, in general, be split 

in ■ 

in a gauge-invariant way into a spin and orbital part. However, the projection of the gluon spin 
onto its direction of motion i.e. its helicity is gauge invariant and is measured in deep inelastic 
scattering on nucleons. The Ji sum rule, relating the quark angular momentum to generalized 
parton distributions, though not based on the canonical operators, is shown to be correct, if 
interpreted with due care. 

I also draw attention to several interesting aspects of QED and QCD, which, to the best of my 
knowledge, are not commented upon in the standard textbooks on Field Theory. 

PACS numbers: 11.15.-q, 12.20.-m, 12.38.Aw, 12.38.Bx, 12.38.-t, 14.20.Dh 
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I. INTRODUCTION 



A major controversy has arisen in QCD as to how to split the total angular momentum 
into separate quark and gluon components (throughout this paper "quark" will mean a sum 
over all flavours of quarks and antiquarks). The idea of identifying separate quark and 
gluon angular momentum operators is attractive, since these operators may be measurable 
in certain physical processes and there may be sum rules relating the spin of a nucleon 
to the angular momentum carried by its constituents. The operators for total momentum 
and total angular momentum, obtained via Noether's theorem from the QCD Lagarangian, 
consist of separate terms which seem to represent a natural division into quark and gluon 
pieces. However, Ji, in particular, [l| has argued that such terms are not individually gauge 
invariant and has advocated use of the Bellinfante version of these operators, which has the 
nice property that they are gauge invariant and can be measured in Deeply- virtual Compton 
Scattering reactions 2]. But Ji's quark angular momentum operator contains both quark 
fields and the gluon vector potential, so is not obviously to be interpreted as the physical 
quark angular momentum. Indeed, a major debate has arisen as to whether it is correct 
to identify this operator as the quark angular momentum, and Chen, Lu, Sun, Wang and 
Goldman [3] and Wakamatsu J4] have proposed quite different identifications, leading to very 
different statements as to what fractions of momentum and angular momentum the quarks 
and gluons carry in the asymptotic limit Q 2 — > 00. In Ji's Bellinfante approach no attempt 
is made to split the gluon angular momentum into a spin part and an orbital part, in accord 
with the long held belief that such a splitting cannot be done in a gauge-invariant way. But 
both Cben at al S and WaWsn Q Cann mu cn m o re , na m e ly , that it is possible to ca rry 
our such a division in a gauge-invariant way and that even in QED the traditional, decades- 
old textbook method of identifying electron and photon angular momentum is incorrect! ( 
For access to the papers in the controversy see ref. [5].) 

The paper of Wakamtsu (4] explains very clearly how the differences between the various 
approaches arise. In QED one splits the photon vector potential into two parts 



corresponding exactly to what is usually called the transverse A± and longitudinal A\\ parts 
respectively, with 

V . Aphya = and V x A pure = (2) 
Under a gauge transformation A p h ys is invariant, whereas 

Ap Ure {x) y Ap Ure {x) + VA(a;) (3) 

In QCD, analogously, one splits 

= 4 M -I- A^ (A) 

a phys,a ' pure,a \ J 

where A^ urea transforms like itself under gauge transformations, but is a pure gauge in 
the sense that it gives rise to no non-zero fields i.e. G^ re = 0, while A^ hys a transforms 
covariantly i.e. like G^ v itself. 

Wakamatsu shows that the difference between the various versions lies in the freedom to 
insert a particular term 

V = g J d 3 x^(x)(x x A a phys )tlJ) m {x) (5) 

either into the quark orbital angular momentum or into the gluon angular momentum, 
yielding, he claims, two possibilities. But, in fact, if there is no other criterion to indicate 
which is the correct choice, there is actually an infinite number of possibilities i.e. one could 
insert aV into the quark orbital term and (1 — a)V into the gluon term. 
In a later paper 6] Wakamatsu attempted to reformulate his approach in a manifestly 
covariant form and to relate his spin and orbital terms to the polarized parton densities which 
are measured in polarized deep inelastic scattering. Unfortunately many of the equations in 
this paper are incorrect as a result of treating a non-forward matrix element like 
(p+A/2; S | M' luX | p— A/2; S ) as transforming like a tensor, and forgetting that the physical 
requirement on the covariant polarization vector, namely S ■ (p±A/2) = implies S ■ A = 0. 
It should be stressed that the existence of these errors is not controversial. The same 
errors occur in the Jaffe-Manohar paper fl aad have been graciously acknowledged by those 
authorsj. Umfortunately then, it is very difficult to decide which claims in the Wakamatsu 
paper are justified. 



1 private communication from Professor Jaffe to T. L.Trueman and the author 
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We shall argue that none of these prescriptions is generally correct or physically plausible, 
but we shall see that the Bellinfante version works in certain specific situations. There are 
three main problems: 

1) In all these papers much emphasis is placed on the issue of using gauge-invariant 
operators. We shall show that this emphasis is misplaced and that the gauge invariance of 
the operators is not an important criterion. In particular we suggest that neither Ji's, Chen 
et al's nor Wakamtsu's identification is physically correct. We shall first show below that in 
any theory which is invariant under gauge transformations, even the total momentum and 
angular momentum operators cannot be gauge invariant. Of course this does not mean that 
the momentum and angular momentum cannot be measured. Because what one measures — 
and this is the key point — are not operators but matrix elements of operators, and if care is 
exercised in defining the physical states of the theory (respecting any subsidiary conditions, 
which is crucial in a gauge theory) then these matrix elements turn out to be gauge invariant. 
This is the basis for our suggestion that the emphasis on utilizing gauge invariant operators 
is misleading. Then we shall discuss what happens if one insists on using gauge invariant 
operators and demonstrate that they do not, in general, have the physical meaning expected 
of them. 

2) In all the above papers the treatment is essentially classical and use is made of the 
classical equations of motion. This totally ignores the highly non-trivial complications in- 
volved in quantizing a gauge theory and the fact that some classical equations cannot be 
maintained at the operator level. For example in QED, when one writes for the photon vec- 
tor potential the symbol A^x), it creates the expectation that it transforms like a 4- vector 
under Lorentz transformations. Yet to agree with the Maxwell equations A^(x) has to sat- 
isfy a subsidiary condition. In classical electrodynamics one chooses the beautiful covariant 
Lorenz condition d^A^x) = 0, which indeed permits A^x) to transform as a 4-vector. It is 
well known, however, that one cannot impose such a subsidiary condition on the operators 
Apix) in QED, since it contradicts the usual canonical equal-time commutation relations of 
the quantized theory. There are many approaches to the quantization of electrodynamics in 
which a non-covariant subsidiary condition is imposed (for a concise summary see Section 
21.2 of jsj]). A popular choice is the Coulomb gauge condition V . A = (see, for example, 
Section 13.5 of |9j) . If this gauge condition is to hold in any reference frame then clearly 
Ap(x) cannot behave as a 4-vector, but — and this is the crucial point — this does not spoil the 
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Lorentz invariance of the theory, since the matrix elements corresponding to any measurable 



physical quantity do transform correctly 10( . Thus, firstly, we suggest that it is unnecessary 
to insist that A^(x) transforms as a 4- vector, and secondly, but more importantly, if, as Ji 
does, one does insist that one's vector potential is a genuine 4-vector, then one has to deal 
with a covariantly quantized theory, in which case the expressions given in the Ji, Chen 
et al and Wakamatsu papers, for the linear and angular momenta, are incomplete. The 



covariant quantization of QED is a non-trivial task 
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13] involving the introduction of a 



scalar gauge-fixing field B(x). Covariant QCD is even more complicated, both in instant 



form 14j and light-front form [15] . involving both a gauge-fixing field and Faddeev- Popov 
ghosts fields. In both QED and QCD the expressions for the linear and angular momentum 
should include terms involving all these fields. 

3) The key issue of splitting the total momentum and angular momentum into a quark and 
gluon contribution is not adequately analyzed. There are two rather separate aspects. There 
is the age-old question of splitting the angular momentum of a gauge particle into a spin 
part and an orbital part. We shall discuss this in Section rVHI But there is a more general 
question of how, in any theory with interacting fields, say <Pe{x) and </>f(x), one can split 
the total momentum (and angular momentum) into pieces interpretable as the contributions 
of the quanta E and F. In all the above papers, having invented some strategy for defining 
the operators Pe and Pp, one writes, for the total momentum 



E 



(6) 



and then interprets the nucleon expectation values of theses operators as a measure of 
the contribution of E and F respectively to the momentum of the nucleon. But this is 
potentially misleading, because the interacting particles constantly exchange momentum, 
and the correct way to express Eq. (J6j) is 

P = P E (t) + P F (t) (7) 

to reflect the fact that while the total momentum is conserved, the individual momenta 
are not. Thus it requires some analysis to explain why it is meaningful to interpret, e.g. 
(nucleon | Pe{^) I nucleon) as a fixed number measuring the contribution of E to the mo- 
mentum of the nucleon. The correct way to extract a measure of the separate contributions 
is to remember, as stressed by Jaffe and Manohar [7( in the QCD case, that constituent 



quark models and parton models of the nucleoli are canonical Fock-space models. Thus 
the physical nucleon states of the theory are taken to be superpositions of quark and gluon 
Fock states. Similarly, in QED, atomic states are regarded as superpositions of electron and 
photon Fock states. How this affects the extraction of the separate momentum and angular 
momentum contributions is spelled out in Section IVII 

We shall argue that the various prescriptions given by Chen et.al. and Wakamatsu are some- 
what ad hoc and that what is missing is a compelling criterion for identifying a particular 
operator as the momentum operator or as the angular momentum operator. The natural 
definition of the total momentum operator is as the generator of translations and of the total 
angular momentum operator as the generator of rotations, but when the system consists of 
different interacting quanta some modification is unavoidable. We suggest that the minimal 
requirement for this identification is the following: 

Definition: Suppose we have a system consisting of interacting fields 0_e(x) and 4>f(x). 
Then the momentum operator P 3 E {t) for, say, particles E should, at equal times, satisfy 

t{P> E (t)^ E (t,x)] = d^ E (t,x). (8) 
Analogously, the angular momentum operator M E (t) should, at equal times, satisfy 

i[M${t) , 0f (t, x)} = - tf&)tf{t, x) + (S^) r s 0f (t, x) (9) 

where r and s are spinor or Lorentz labels and (E 1 -?)/ is the relevant spin operator. The 
need for the requirement u at equal times" is explained in detail in Section IVT1 
Demanding that these conditions be satisfied leads to the conclusion that the canonical 
expressions for the momentum and angular momentum operators are the correct and phys- 
ically meaningful ones. It is then an inescapable fact that the photon and gluon angular 
momentum operators cannot, in general, be split in a gauge invariant way into a spin and 
orbital part. However, as discussed in Section PVIIl the projection of the photon and gluon 
spin onto their direction of motion i.e. their helicity, is gauge invariant and is measured in 
deep inelastic scattering on atoms or nucleons respectively. 

It should be noted that Ji's expressions for the components of the quark and gluon momen- 
tum and angular momentum vectors, which are the Bellinfante versions, do not conform to 
the above definition and thus should not be considered as measuring all the components of 

6 



the physical quark and gluon momentum and angular momentum vectors, though it turns 
out that they give the correct results for the Z-components, P z and J z , for a nucleon moving 
in the Z direction. In particular the quark orbital angular momentum defined by Ji as 
the difference between his quark total angular momentum, as measured in Deeply-virtual 
Compton Scattering, and the quark spin, as measured in Polarized DIS, is in agreement 
with our definition, as long as it is appreciated that this refers only to the components along 
the direction of motion of the nucleon. 

The difficulty in defining separate quark and gluon angular momenta in QCD has its ana- 
logue in QED, in the problem of defining separate electron and photon angular momenta. 
However, the situation is not completely analogous in the two cases, because the straight- 
forward gauge invariance of QED is replaced by the rather different BRST [if]] invariance 
of QCD. For this reason we shall discuss the two cases separately. 

Most of the problems which beset the definition of separate quark and gluon angular 
momenta actually already occur at the level of the linear momentum. Since this is a much 
simpler object to deal with, we shall mainly illustrate the problematic issues through an 
analysis of the linear momentum operator. 



II. OBSERVABLES IN GAUGE THEORIES 

As mentioned above we think there has been too much emphasis on the need to use gauge 
invariant operators to represent any dynamical quantity which can be measured i.e which 
is an observable. In this section we shall show that, in fact, in gauge theories the concept 
of an observable is very subtle and is rather different in QED and QCD, and we shall give 
the precise conditions that an observable operator must satisfy. Our discussion follows the 



approach o 
Wightman 



Kugo and Ojima [17j], which, in turn, follows the treatment of Strocchi and 



In the covariant quantization of a gauge theory it is unavoidable that one has to set up 
the theory in a vector space with an indefinite metric i.e. one in which the "length" or norm 
of a vector can be negative. From this one constructs a subspace, the physical vector space 
Vptys, in which scalar products are positive semi-definite, and finally the positive definite 
Hilbert quotient space Uphys = Vphys/Vb, where Vo is the subspace of V p h ys consisting of 



zero-norm vectorgj. How the states forming V p h ys are defined, depends on the formulation 
of the theory. In covariantly quantized QED they are defined by B^^x)]^) = 0, where 
B{x) is the gauge fixing field. In covariantly quantized QCD one has = , where 

Qb is the generator of BRST transformations. 

Note that while V p h ys is labelled "phys" , the states which correspond to the actual physical 
particles belong to Uphys i-e. the zero-norm states in V p h ys are not truly physical. We assume, 
as usual, that the physical states form a complete set in Uphys- 

Let |$) be a state in V p h ys i-e. |$) G V p h ys and let \x) be a zero-norm state i.e. \x) G Vo- 
Then it can be shown that 

($| X ) = for V|$)GV^ S , V|x)GV (10) 

i.e. 

Vo-L V phys . (11) 

Let O be a physical quantity and let O be the hermitian operator representing it. It can 
be shown that a necessary condition for O to be an observable is 

+ + x ) = ($|6|$) for v\$)ev phys , V|x)GV (12) 

Equivalently, via Eq. (II ip . an observable operator must satisfy 

0\$)e V ph ys for V|$)G Vphys (13) 

The essential point of this argument, as we shall see later, is that the condition Eq. ( TT3|) 
does not necessarily require an observable operator to be gauge invariant in the operator 
sense i.e. to commute with the generator of gauge transformations. And we shall see that 
the situation differs somewhat between covariantly quantized QED and QCD. 

2 Strictly speaking 7i p hys should be defined as the completed quotient space, but this is irrelevant for our 
discussion. 
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III. THE MOMENTUM OPERATOR IN GAUGE-INVARIANT THEORIES 



If the theory is invariant under translations in space-time, then Noether's theorem allows 
the construction, from the classical Lagrangian, of what is usually referred to as the canonical 
energy-momentum tensor density t^ n (x). This is a conserved density 

d,t c Z(x) = (14) 
but is generally not symmetric under //•<-)■ v. 

The canonical total linear momentum operator P-! an is the space integral 

PL = j d 3 x t°J n (x) (15) 
and, crucially, is independent of time as a consequence of Eq. f JT4|) . 

A. The canonical momentum operator as generator of translations 

In the classical theory P^ an thus constructed is the generator of spatial translations. In 
the quantum theory one has to check that the operator version of P-! an satisfies the correct 
commutation relations with all the fields i.e. for any field (j)(x) 

i [P^ cj>{x)] = &4>(x) (16) 

It is important to realize that in an interacting field theory an arbitrary commutation 
relation between the fields cannot be calculated unless one can completely solve the theory — 
an impossible task in all relevant physical theories. On the other hand the Equal Time 
Commutators (ETC) are fixed as part of the process of quantizing the theory. Hence the 
only reason it is possible to check an equation like (JTB]) is because P^ an is independent of 
time and so the time variable in the fields occurring in it can be chosen to coincide with the 
time variable in <j)(x) = <p(t, x). This consideration will play a crucial role when we come to 
discuss how to divide the total momentum into contributions from the different fields in the 
theory. 

An important issue in comparing the treatment of linear and angular momentum is the 
concept of a local operator. An operator 0{x) is local if, obviously, it is defined at one 
space-time point x, but also it must satisfy the law of translation 
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0(t,x + a) = e iPL " a > 0(t, x) e' iPLn a = . (17) 

Note that an operator of the form M(x) = xO(x), such as occurs in the expression for the 
angular momentum, is not a local operator. (It is trivial to see that if M(x) satisfies Eq. (TT7T) 
then M(x) = for all x.) In a careful discussion of the properties of angular momentum, 
operators of this type have been called compound operators 191 ] . 



B. The Bellinfante energy momentum operator tensor density 

As mentioned the canonical t^ n (x) is generally not symmetric under interchange of \x and 
v. It is also not gauge invariant. It is possible to construct from t^ n (x) and the Lagrangian, 
the conserved Bellinfante density t^(x), which is symmetric and, which is, in some cases, 
as will be discussed below, gauge invariant. It differs from t^ n [x) by a divergence term of 
the following form: 

O) = &) + \d p [H"» - H"» - H»n (18) 

where the only relevant property of H ppu for the present discussion is that it is antisym- 
metric under \i <H- v 



and that it is a local operator. 
It follows that P d hel defined by 

PL = J d 3 x C(x) (20) 

differs from P^ an by the integral of a spatial divergence, and it is usually stated that 
since the fields must vanish at infinity, such a contribution can be neglected, leading to the 
equality 



^bel ~ Plan- (21) 

Now for a classical c-number field it is meaningful to argue that the field vanishes at 
infinity and that Eq. Q2ip holds as a numerical equality. It is much less obvious what this 

10 



means for a quantum operator. The correct way to tell whether a divergence term can be 
neglected is to check what its role is in the relevant physical matrix elements involving the 
operator. In the case of Eq. (1211) one can readily check that the matrix elements between 
any normalizable physical states, |^) and |$) are the same^ i.e. 



<*|i&|¥) = <*l ( 22 ) 

However, the operators cannot be identical, because one, for example, may be gauge 
invariant and the other not, so that the equality would be contradicted upon performing a 
gauge transformation. On the other hand the operators are essentially equivalent, and they 
generate the same transformations on the fields. We shall indicate the relationship as 



Pbel — ^can- (23) 

It should be noted that it would be impossible to construct a consistent theory if it 
were not permissible, in certain case, to ignore the spatial integral of the divergence of a 
local operator. For example we could not even establish the obvious requirement that the 
momentum operator commutes with itself! For one has, (no sum over j) 

i[P j , P j ] = J d 3 x i[P j , t 0j (x)} = J d 3 x dH 0j (x) (24) 

and this vanishes only if the divergence integral can be ignored. 

For compound operators like the angular momentum it is a much more difficult task to 
show the equivalence of the total angular momentum generators M l J an and M^ z , constructed 
from the canonical and Bellinfante P ca n,bei respectively, and care has to be exercised to always 
use normalizable states. This has been done by Shore and White 20]. 



IV. QUANTUM ELECTRODYNAMICS 

We shall study the questions of gauge invariance and Lorentz covariance first in the 
simpler context of QED. 

3 This is not true for all operators which differ by a divergence term. Singularities can affect the result. 
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A. The non-gauge invariance of the QED momentum and angular momentum 
operators 

We remarked in the Introduction that in trying to define separate quark and gluon angular 
momentum operators too much emphasis was being placed on the use of gauge invariant 
operators by Ji, Chen et al and Wakamatsu. 

In support of this point of view we shall now prove that in any theory which is invariant 
under a local c-number gauge transformation, even the total momentum and angular mo- 
mentum operators cannot be gauge invariant. As discussed in Section [TT] this does not mean 
that the momentum and angular momentum are not observables i.e. cannot be measured. 
Because what one measures are not operators but matrix elements of operators, and if care is 
exercised in defining the physical states of the theory (respecting any subsidiary conditions) 
then these matrix elements turn out to be gauge invariant. 

Theorem 1: Consider a theory which is invariant under local c-number gauge transfor- 
mations. Let P^ be the total momentum operators, defined as the generators of space-time 
translations, and let M 1J be the total angular momentum operators, defined as the genera- 
tors of rotations. Then P M and M u cannot be gauge invariant operators. 
Proof: For simplicity we consider QED and give the proof just for the momentum oper- 
ators. The case of angular momentum is a straightforward generalization. Note that it is 
irrelevant for the proof whether we use the canonical or Bellinfante versions. 
The theory is invariant under the infinitesmal gauge transformation 



A"(x) ->■ A"(x) + d»A(x) 



(25) 



where A(x) is a c-number field satisfying □A(rr) = and vanishing at infinity. 
Now gauge transformations are canonical transformations 2l|. Let F be the generator 
of gauge transformations, so that 



i[F,A^(x)) — d' x A(x) (26) 

and consider the Jacobi identity 



[F, [P», A"}] + [A", [F, P»)) + [P», [A", F}) = 
12 



(27) 



Now [P M , [A^P]] = since by Eq. (|26|) [A^F] is a c-number and thus commutes with 
P", so that 



[[F, P% A"} = [F, [P^A V \] (28) 
Moreover since P^ are the generators of translations 

i[P»,A u ] = d ft A'' (29) 
Thus the RHS of Eq. (|2B|) becomes 

[F, [P", A"]] = A v {x)\ = -d"d v K{x) ^ (30) 

and hence from Eq. ([28"]) 

[[P,P1,A1^0 (31) 

implying that 

[F, P"} ^ (32) 

so that P M is not gauge invariant. 



B. The momentum and angular momentum in QED are observables 

We shall now demonstrate that this lack of gauge invariance is of no physical significance. 
We shall take as an example covariantly quantized QED and show that the matrix element 
of Pl an between any physical states, is unaffected by gauge changes in the operator. 

As far as we are aware the most general covariantly quantized version of QED is given by 



the Lautrup-Nakanishi Lagrangian density III, ll2|, which is a combination of the Classical 



Lagrangian (Clas) and a Gauge Fixing part (Gf) 



£ — £das + &Gf (33) 

where 
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— m + e JL)tp + h.c. 



(34) 



and 



£ Gf = B(x)d tl A»(x) + -B 2 (x) 



(35) 



where B(x) is the gauge-fixing fielcj^ and the parameter a determines the structure of 
the photon propagator and is irrelevant for the present discussion^. The theory is invariant 
under the usual c-number infinitesmal gauge transformation 



— >■ + d^A(x) ip -»■ V + ^eA^ (36) 

while -B(x) is taken to be unaffected by gauge transformations. 

A straightforward calculation gives for the conserved generator of infinitesmal gauge 
transformations 



F = - J d 3 x [e^7°#(x) + F° j djA(x) - B(x)d A(x)} 
which,via the equations of motion, can be transformed to 

F = J <fx [(d B)A - Bd A + <9,-(F oi A)]. 
Now the physical states |<E>) of the theory are defined to satisfy 



(37) 



(38) 



= 



where 



(39) 



fl(i)=Bl +) (i)+fl ( "l(i) (40) 
with B^'fi) the positive/negative frequency parts of B(x). 

With this definition of the physical states, an operator O is an observable, if, according 
to Eq. (fl2l) . 0\ $ ) is itself a physical state i.e. if 

fi (+) (x)(6|$» = 0. (41) 
4 Because of its similarity with the QCD case, we use the notation of Nakanishi. Note that Lautrup's 



A(a:) = -B(x). 

5 The case a = 1 corresponds to the Gupta-Bleuler approach (see e.g. 



21]) based on the Fermi Lagrangian. 
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This is equivalent to the condition 

[B {+ \x), 6]|$) = (42) 

since 

[B M (x), 6]\$) = B {+ \x)0\<5>) -OP (+) (x)|$) = B M (x)d\<£>) 

= iff Eq. (JUD holds. (43) 

Since, via Eqs. (ESJ and (USD 

[B (+) (x),F] |$) = ^'S (+) (x) |$) = (44) 

we see that P J is an observable, so that its eigenstates are physical states. 

We shall now consider the gauge invariance of its matrix elements. In doing so it should 
be noted that B { -\x) = (x), so that )(.x) = 0, and thus for arbitrary physical 

states 

(&\B(x)\&) = 0. (45) 

Theorem 2 Any physical matrix element of the momentum operator P- 7 is invariant 
under gauge transformations. 

Proof Consider the general physical matrix element 

($'|pj|$) = J d 3 pd 3 p'(f)'*(p')4>(p)(p'\P j \p) (46) 

The change induced in (p'l-P-'lp) by the gauge transformation is given by (p'\i[F, Pi}\p). 
Focus initially on the effect of the first two terms (call them f\ 2 ) in the integrand on the 
RHS of Eq. flSE]). 

(p'\i[fi2,P j ]\p) = (P-P'Y (p'lfM = (47) 

because of Eq. (T45|) and the fact that A is a c-number. 

The change induced by the third, divergence term (call it Z^) in the integrand on the RHS 
of Eq. (I3"81) . after some algebra, and using translation invariance Eq. ffT7|) . can be written 



(pWs,P J ]\p) = {p'-py{{p-p'f{p'\A k m P ) (48) 
- (p-p'f <jpf\A\0)\p)}d k [A(x)e*>^*] (49) 
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and this vanishes after the spatial integration because A(x) vanishes at infinity. 
Hence ($'|P J '|$) is indeed invariant under gauge transformations. 

Corollary Physical matrix elements of the total angular momentum operator J k are 
gauge invariant. 

The total angular momentum is given by 

J k = \e klm M lm = l -e klm j d 3 xM 0lm (x) (50) 

where A4 0lm (x) is the angular momentum tensor density. The simplest way to show the 
gauge invariance of the physical matrix elements in this case is to reinterpret the gauge 
change i[F, M lm ] as —i[M lm , F] and to study the effect of the rotations on F. For this one 
needs the following results: 

i[M lm , AP(x)\ = {x l d m - x m d l )A /3 (x) + g l/3 A m (x) - g m/3 A l (x) (51) 

and, since B(x) is a scalar field, 

i[M lm , B{x)\ = {x l d rn - x m d l )B(x) (52) 

Application of these to F yields terms which either vanish directly as a result of the sub- 
sidiary condition Eq. (T4"5]) or divergence terms which can be shown to vanish since A(x) 
vanishes at infinity. 

The fact that even the total momentum and angular momentum are not gauge invariant, 
but that their physical matrix elements are, suggests that to insist on gauge-invariant op- 
erators for the momentum and angular momentum operators of the individual fields of the 
theory is unnecessary. 

C. Relativistic covariance in QED 

In the debate with Chen et al, Ji rightly argues that their photon vector potential does 
not transform as a 4- vector under Lorentz transformations, and implies that in his treatment 
his A^{x) transforms as a true 4- vector, and that this is an essential property. But if this is 
the case then Ji's expressions for momentum and angular momentum are incomplete. The 
point is that the gauge-fixing field B{x) introduced above in the covariant quantization of 
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QED also appears in the expressions for the momentum and angular momentum. One finds 
for the conserved canonical energy momentum tensor density, 



C = C + C(G/) (53) 

where 

= \ 1> ~ FrtdTAf, - g^Ccus (54) 

where ^ v = — *d u , and 

tZn(Gf ) = Bd»A» - g^C Gf . (55) 
For the conserved Bellinfante density one finds, 

C = C + C(G/) (56) 
where 9^, which is referred to as the classical energy momentum tensor density, is 

C = \ tW^" + Y^") 1> - Frt^ - gTCoas (57) 

^ y ^ y 

where D v — d v - 2ieA u , and 

€i(Gf) = -(d»B) A" - {d v B) A" - g^C Gf (58) 

The conservation of an energy momentum tensor depends on the equations of motion, 
which are a consequence of the Lagrangian. Thus t^ el is conserved, but 9^ is not, when the 
Lagrangian is Ccias + C-Gf- On the other hand 9^ would be conserved if the Lagrangian 
were C C i as - 

Now Ji and Chen at al utilize 9^ v el and treat it as if it were conserved i.e. they take 
the momentum operator based on it to be independent of time (equivalently: to remain 
unrenormalized) , which implies that the Lagrangian is just Ccias- But it is well known that 
one cannot quantize QED covariantly using C = Hcias- Nonetheless Ji insists that his A M 
transforms covariantly, which is thus, at the operator level, a contradiction. 

We have seen that insisting on covariant quantization leads to a more complicated struc- 
ture for the energy momentum density and analogously for the angular momentum density. 
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This raises what, at first sight, seems to be a worrying issue concerning several papers in the 
literature, e.g. Ji [l, 2|, Jaffe and Manohar [7], Bakker, Leader and Trueman (BLT) [19] and 
Wakamatsu |4|, |6j], where the general structure of the matrix elements of 9^ (or its QCD 
analogue) is derived under the assumption that 8^ is a genuine conserved tensor. However 
the situation is saved by the following: for physical matrix elements, for both the canonical 
and Bellinfante versions, 



($'| t^ u (Gf) |$) = 0. (59) 

This follows from Eqs. fl5"5"l |5"S|) and (1551) when a complete set of physical states is inserted 
between the operators appearing in t^ u {Gf) and use is made of Eq. (H5l) . Hence 

m <v») i*> = <*i i*> = o. (60) 

Similar arguments show that 0^ n {x) , which just corresponds to the canonical version of 
0^(x), may also be treated as a conserved density inside physical matrix elements. Thus 
this aspect of the analysis in the above papers is, in fact, consistent. 

In summary covariant quantization of QED complicates some aspects and there is no 
compelling reason to insist on it. Indeed, as explained in the Introduction, the non-covariant 
Coulomb gauge leads to a perfectly good Lorentz invariant theory. However, if one prefers 
to work with a covariantly quantized theory then, in so far as its physical matrix elements 
are concerned, #^(x) and 0^ n (x) may be treated as conserved tensor operators. 



V. QUANTUM CHROMODYNAMICS 

The situation in QCD is somewhat different. The infinitesmal gauge transformations on 
the gluon vector potential and on the quark fields, under which the pure quark-gluon La- 
grangian C qG ( the QCD analogue of the QED Ccias), 

Ac = -\g%G? + \^%mi{$ - V) - 2< m A a ]^ m (61) 
is invariant, are determined by eight scalar c-number fields 8 a (x), 



5A« = d,6 a (x)-gf abc AZ(x)9 c (x) 
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(62) 



= -igtt m a (x)Tp m {x) (63) 
where a,b,c — 1, 2. ..8 and I, m = 1, 2, 3 are colour labels, and where our sign convention 

is 



G% = d,A a v - d v Al - gfabcA^Al. (64) 

However, in order to quantize the theory covariantly one has to introduce both a gauge- 
fixing field B(x) and Fadeev-Popov anti-commuting fermionic ghost fields c(x), c(x). The 
Kugo-Ojima Lagrangian [14| for the covariantly quantized theory is then 

£ = C-qG + £-Gf+Gh (65) 

where 

C Gf+Gh = -i(d»c a )Dfc b - (d»B a )A* + ^B a B a (66) 

which is no longer invariant under the original infinitesmal gauge transformations 
Eqs. ([621 [63]). 

One can again show that the momentum operators P can , Pbei are not gauge invariant, 
but this is now irrelevant, given that the Lagrangian itself does not possess this invariance. 
Instead the theory is invariant under the BRST transformations jl^] 



6A« = BDfc\x) 

= -i6gt? m c a (x)ij m (x) 
6c a (x) = 9(g/2)f abc c b (x)c c (x) 
Sc a = i6B a (x) 

5B(x) = 0. (67) 

where 9 is a constant operator which commutes with bosonic fields and anti-commutes 
with fermionic fields. 

The BRST transformation is generated by 6Qb i-e. for any of the above fields 4> 
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i[6Qb,<P] =H 
where the conserved, hermitian charge Qb is given by 



(68) 



Q B = J d 3 x[B a ^ c a - gB a f abc A b Q c c - i{g/2)(d c a )f abc c b c c ]. (69) 
There is also a conserved charge 

Q c = j d 3 x[c a f c a - gc a f abc A b c c ] (70) 
which "measures" the ghost number 

t[Q c ,<f ) } = N<j ) (71) 

where N — 1 for = c, —1 for = c and for all other fields. 
The physical states are defined by the subsidiary conditions 

Qb\^} = (72) 

Qc\*) = (73) 

A. The momentum and angular momentum operators in covariant QCD 

The proof of an analogue of Theorem 1 for BRST transformations does not work, 
because the BRST 5A^ is an operator, not a c-number. Consequently, use of the Jacobi 
identity Eq. (1271) . with F replaced by Qb, does not imply that P can or P be i are non-BRST 
invariant. 

Analogously to condition Eq. (142 p . in order to be observable the momentum operator in 
QCD must satisfy 

[Q B ,P j m = and [Q C ,P J M = 0. (74) 

The latter, as will be seen presently, follows from the fact that the ghost number of P J is 
zero. The former is usually stated to hold because Qb is a translationally invariant scalar. 
This is correct, but is not quite as trivial as it seems, for if we write 

Qb = J d 3 xQ B (t,x) (75) 
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then translational invariance requires 

e ipia iQ B e- ipSa * = j d 3 xe ip3a >Q B (t,x)e- iP ' a > = J d 3 xQ B (t,x + a) = Q B . (76) 

The last step holds only if the integral in invariant under the change of variables x — > y = 
x + a, which is in accord with our being able to ignore the integral of a divergence. 
One finds for the canonical energy momentum tensor density, 

Cn = + Cn(G7 + Gh) (77) 

where 

C n (qG) = \frf*$ v 1>i - Gfd»A% - g^C qG (78) 

and where 

W + Gh) = -A£PB a - i{d v c a ){D» b c b ) - gTCaf+GH - i{d»c a ){d v c a ). (79) 
The Bellinfante version is 

d = €i(QG)+tZ(Gf + Gh) (80) 

where 

tZ(QG) = + (a* v)\ ~ GfCfy - g^C qG (81) 

^ y 

is BRST invariant, i.e. commutes with Qb- Here D v is a matrix in colour space 

*d»{z) = 5 lm [S» - + 2igA:(z)tt m . (82) 
The gauge-fixing and ghost terms are given by 
tZ(Gf + Gh) = -{A*PB a + A^B a ) - z[(d»c a )D: b c b + {d v c a )D^c b ] - g^C Gf+Gh . (83) 
This can be rewritten as an anti-commutator with Qb 

C el (Gf + Gh) = -{Q B , {{d»c a )A v a + (d u CaW a + g^c a B a - (VcJAfl)}. (84) 



It follows that t^(Gf + Gh) is BRST invariant (because is nilpotent i.e. Q\ = 0) 
and does not contribute to physical matrix elements i.e. 

(#ICI*> = <#IC(?G9I*>- (85) 
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Thus the entire t^(x) commutes with Qb and is therefore a local observable. 

The situation with t^ n (x) is somewhat different. It does not commute with Q B , so is 



not itself an observable, but, contrary to the statement in [20[, t^ n (Gf + Gh) does not 
contribute to physical matrix elements. This can be seen as follows. The first three terms 
in Eq. (1791) can be written as an anti-commutator with Qb, so, as argued above, do not 
contribute to physical matrix elements. For the last term we have, by completeness, 

-Z<$'| {d»CaWc a ) |$) = ( d ^a) |*><tf| (d V C a ) |$>. (86) 

This is zero because, via Eq. (17T1) . c a (x) = i [Q c , c a (x)], so that 



(tf | d u c a {x) |$> = i9"<*| [Q c , c a (x)] |$) = (87) 
as a consequence of Eq. (1751) . 

Thus even though the actual canonical density is not BRST invariant, its ghost and 
gauge-fixing terms do not contribute to physical matrix elements. And, as discussed in 
Section II III B, for the space integrated versions, because they differ by a divergence, we 
have, analogous to Eq. (120]) . 



PUQCD) - PUQCD) (88) 

and both are BRST invariant. 

There is thus no compelling reason in QCD for insisting on using the Bellinfante version. 
Analogous statements hold for the angular momentum generators M l J an and M^. 



B. Relativistic covariance in QCD 

We have seen that insisting on covariant quantization forces us to include gauge-fixing 
and ghost fields in the Lagrangian. However, the terms in the canonical and Bellinfante 
versions of the total momentum, which depend on the ghost and gauge-fixing fields, do 
not contribute to physical matrix elements. Thus if we consider the expectation value of 
the total momentum operator for a nucleon in a state of definite momentum | p ) then, 
irrespective of whether we use P can or P beh there will be no contribution from the ghosts or 

22 



gauge-fixing fields. Moreover, both operators are observables and their matrix elements are 
thus physically measurable quantities. 

However, as mentioned in the Introduction, Ji, Chen et al and Wakamatsu insist on using 
the gauge- invariant Bellinfante tensor, or modifications of it, for the separate electron and 
photon, or quark and gluon, parts of the total momentum and angular momentum tensors. 
We shall argue in the next Section that this has no solid basis, is essentially arbitrary and 
lacks any persuasive physical motivation. 

VI. THE PROBLEM OF DEFINING SEPARATE QUARK AND GLUON MO- 
MENTA 

We come now to the heart of the controversy between Ji, Chen at al and Wakamatsu, 
namely how to define in a sensible way the separate contributions of quarks and gluons 
to the momentum and angular momentum of a nucleon. There are actually two separate 
issues. One, quite general, is how to define the separate momenta for a system of interacting 
particles. The second is more specific to gauge theories and includes the issue of splitting 
the angular momentum of a gauge particle into a spin and orbital part. 

A. Interacting particles: the general problem 

Suppose we have a system of interacting particles E and F and we split the total mo- 
mentum into two pieces 



which we wish to associate with the momentum carried by the individual particles E and 
F respectively. 

As mentioned in the Introduction it is crucial to realize that Eq. ( )89|) . as it stands, is 
totally misleading, and should be written 



to reflect the fact that the particles exchange momentum as a result of their interaction. 



P j = P J E + P J F 



(89) 



Pi = P> E (t) + P^t) 



(90) 
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The key question is: what should be the criterion for identifying Pe,f as the momentum 
associated with particles E, F respectively? 

The seductively obvious answer would be to demand that 

i[Pl ( j ) E {x)} = d^ E {x) (91) 

and similarly for F, but there is no way we can check this, since P E (t) depends on t and, 
without solving the entire theory, we are only able to compute equal time commutators . 

We suggest, therefore, that the minimal requirement for identifying an operator P 3 E with 
the momentum carried by E , is to demand that at equal times the analogue of Eq. (I9~T|) 
holds i.e. 

i[Pi(t),<f> E (t,x)} = d^<j> E (t,x). (92) 

Analogously, for an angular momentum operator we suggest the minimal requirement 
that 

i[M${t) , <j> E (t, x)\ = (x^ - x^d^fit, x) + (£^) r >f (*, x) (93) 

where r and s are spinor or Lorentz labels and (S 1 - 7 )^ is the relevant spin operator. 

Now we explained in Section IHIB that for the total momentum there is no essential 
difference between P can and Pb e i, since their integrands differ by the spatial divergence of a 
local operator. However, if we split P can into P can ,E + P can ,F and P hei into P he i,E + Pbd,F, 
then the integrands of P can ,E and Pb e i,E do not differ by a spatial divergence, and hence 
Pcan,E and Pbei,E do not generate the same transformation on (p E (x), and similarly for F. 

As an example consider QED. From Eqs. ( l53l 1541) and Eqs. ( |56l l57l l58l) we would identify 

^(electron) = |fy°^> (94) 

and 

tgj (electron) = % - ^( 7 ^' + 7 ^°) ^ (95) 
and these do not differ by a spatial divergence. 

It should be noted that the difference between various definitions of the momentum 
operators is not just a question of principle. In QCD the asymptotic (Q 2 — > oo) limit 
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of the longitudinal momentum carried by quarks in a nucleon, with the Ji definition is 
P(quarks)jj ~ 50% whereas with the Chen et al version P(quarks)c/ien ~ 80%, for the 
number of flavours rif = 5. 

Since, by construction, P can ,E and P ca n,F generate the correct transformations on <Pe{x) 
and respectively, we conclude that with the above minimal requirement we are forced 

to associate the momentum and angular momentum of E and F with the canonical version 
of the relevant operators. This disagrees with Ji, Chen et al and Wakamatsu, but agrees 
with Jaffe and Manohar 7|. 

Nonetheless, exceptionally, for the fraction of the Z-component of the longitudinal momen- 
tum and angular momentum carried by the quarks in a nucleon moving in the Z direction, 
the distinction between Bellinfante and canonical versions is not crucial, since it turns out 
that P 2 (quarks)ji = P 2 (quarks) be ^ = P 2 (quarks) can and J 2 (quarks)ji = J z (quarks) hei = 
J 2 (quarks) can , as will be discussed in Section [VI El 

Now, as Jaffe and Manohar [7] have emphasized in the QCD case, constituent quark 
models and parton models of the nucleon are canonical Fock-space models. Thus the physical 
states of the theory are taken to be superpositions of Fock states, formed from the vacuum 
by the quark and gluon "in-field" creation operators. Similarly, in QED, atomic states are 
regarded as superpositions of Fock states, formed from the vacuum by the electron and 
photon "in-field" creation operators. Loosely speaking, for any field <j)(x) ^ 

(j){x) VZ(j) m (x) (96) 

where Z is a renormalization constant. Also 

[H^bWiMI*)] 1 ^ [PL,can(E),VZ(l>in,Efrx)} (^) 

&<j> E (t, X)^^ VZ&^Eit, x) (98) 

where we have defined 

PL,e(*) ^> Pican(E). (99) 

Note that because the "in" fields obey free field equations, Pf n can (E) is independent of time. 

Now as we have stressed Pb,f(^) are time-dependent operators. However, these operators 
possess a remarkable property. While their general matrix elements are time-dependent, 
6 Strictly speaking such limits of operators should be carried out using normalizable "smearing functions" . 
We shall continue to be a little cavalier in order not to complicate the presentation. 
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there is a sub-class of these, and it is just this class of matrix elements that are of interest 
to us, which are time-independent, namely, their matrix elements between arbitrary states 
of a single particle. To see this for the momentum P E (t) let 

|^) = J d 3 p'ip(p')\p') and |0) = J d 3 p<f)(p)\p) (100) 
be arbitrary states of of a particle of mass m, so that 

pl=p 2 + m 2 and p% = p' 2 + m 2 . (101) 

Then 

<v>|Pi(oi0) = J d 3 p'd 3 pd 3 xr(p > )^p)(p > \t^(x)\p} 

= J d 3 p'd 3 pd 3 x^*(p') <j)(p)e ix - {p '- p) e ^-Po)(p' | tg'(0) \p) 
= (2tt) 3 J d 3 p'd 3 p^*(p')(f)(p) 5 3 (p' -p)e it{p °- p 'o\p'\t o J(0)\p) 

= (2n) 3 J d 3 P r( P )<i>(p)(p\t M\p} (102) 
which is independent of time because p' = po = \/ p 2 + m 2 . 

A similar, though more complicated argument, shows that the single particle matrix 
elements of the angular momentum operators J % E F are also time-independent. 
It follows that e.g. 

(lM^(f)l^>= lim | Hit) \<t>) = (i>\PL(E)\<P) (103) 

t—>— CO 

and analogously for the angular momentum operators. 

Thus we have the important result that the nucleon matrix elements of P 3 E F and J E F can 
be calculated by inserting a Fock expansion for the nucleon state and then evaluating the 
Fock state matrix elements of the "in" field operators P- n (E), P- n (F), J 3 in {E) and J 3 in {F) 
respectively. 

B. Interacting particles in gauge theories: canonical vs "the rest" 

The objection of Ji, Chen et al and Wakamatsu to the use of the canonical operators is 
that they are not gauge invariant. We have suggested that this is not obviously important 
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since the total canonical momentum and angular momentum operators are observables and 
their physical matrix elements are gauge invariant or BRST invariant (Sections IIVB . IVR). 
That argument relied on the fact that an arbitrary physical state can be expressed as a 
superposition of eigenstates of total momentum. 

Now from Eqs. (19611991) and Eq. (|9~2"j) it follows that, for E (and analogously for F) 

i[PL,can(E), <f>in,E(t^)} = & *) ( 104 ) 

which implies that the Fock states, created from the vacuum by the action of the creation 
operators in <f) in>E (x), (f) in)F (x), are eigenstates of P ] in>can {E) and Pj n>can {F) respectively. 
This fact will be used in the next two sections in proving the gauge or BRST invariance of 
the Fock space matrix elements of P 3 n can (E) and P 3 n can (F). An analogous statement holds 
for the angular momentum operators. 



C. QED 

Here particles E and F correspond to electrons and photons and the Fock states may be 

taken as superpositions of states with electrons having definite momentum Pi,p2 p n 

and transverse photons with momenta fei, k 2 k m . It is possible to show that these 

eigenstates of P 3 n can (electron) and P 3 n can (photon) are physical states i.e. 



B in ( x )\Pl,P2 Pn] k 1: k 2 km) = 0. 



(105) 



This follows from the asymptotic limit of the commutation relations given in 11], and the 
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231 ]. which states that the commutators of asymptotic 



Greenberg-Robinson theorem 
fields are c-numbers. 

Since we are only concerned with physical matrix elements of the momentum operators 
we may, as a consequemce of Eq. (11051) ignore the gauge-fixing terms and from now on utilise, 



P 3 

in, can 



(electron) = / d 3 x 



7T Anl° 9 3 Ipi 



and 



P L can (Photon) = / d 3 X 



-F o0 d 3 A 



in, ft 



(106) 



(107) 



The proof that the Fock space matrix elements of these operators are gauge invariant 
requires that the matrix elements of B in (x) vanish between these states. This follows 
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from Eq. (11051) and thus the proof of the gauge invariance of the expectation values of 
P- n can (electron) and P- n can (photon) can be carried through in the same way as was done 
for the total momentum in Section HVl B. 

Note that the simplified versions of the canonical momentum operators above generate 
the correct transformations on ipi n {x) and the spatial components A\ n (x), namely 

i [PL, can (electron), ip in (t, x)} = d J ip in (t, x) i{PL tCan (photon), A^ n {t, x)} = d J A^ n (t, x). 

(108) 

On the other hand one can show that the Bellinfante versions P- n bel (electron) and 
PL & e «(ph°t on ) do not generate the transformations Eq. (11081) . Thus the Bellinfante ver- 
sions do not satisfy our minimal requirement for identifying these operators as representing 
the momentum carried by the electrons and photons respectively. The same is true of the 
Chen et al and Wakamatsu momentum operators. 

The analysis of the angular momentum operators is quite analogous and one concludes 
that the canonical operators are the ones that generate the correct rotations on the fields. 



D. QCD 

Similar results hold for QCD. The states with quarks having definite momentum Pi,P2~ 

p n and transverse gluons having momenta fei, k 2 k m are eigenstates of P- n mn (quark) 

and P- n can (gluon) and are physical states i.e. 



Qb\Pi,P2 



Pn] &1) &2 



km) =0 



(109) 



This follows from the commutation relations for the asymptotic fields given in Section IV in 
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Since we are only concerned with the physical matrix elements of the momentum operators 
we may, as a consequence of the discussion following Eq. ( 18"5|) . ignore the gauge-fixing and 
ghost terms and from now on utilise 



P /n,mn(q Uark ) 



d 3 x 



and 



P 3 ( 

in, can 1 



gluon) = / d 3 x 



^ r in,a u ^in,l3 



(110) 



(111) 
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These commute with Qb and are thus observables. Moreover these simplified versions of 
the canonical momentum operators generate the correct transformations on ip l in (x) and the 
spatial components Af n a (x), namely 

*[^U„(quark), i>\ n t, x)] = &4 n (t, x) i[Pi ca „(gluon), A^Jt, x)] = &A* n Jt, x). 

(112) 

On the other hand one can show that the Bellinfante versions P- n bel (quark) and P( n bd (gluon) 
do not generate the transformations Eq. fll 12j) . Thus the Bellinfante versions do not satisfy 
our minimal requirement for identifying these operators as representing the momentum 
carried by the quarks and gluons. Similar remarks apply to the Chen et al and Wakamatsu 
operators. 

Similarly, one sees that the correct rotations of the fields are generated by the canon- 
ical versions of the angular momentum operators, which suggests that the Ji, Chen et al 
and Wakamatsu operators should not be regarded as representing the angular momentum 
of the quarks and gluons. Nonetheless, the expectation value of the Bellinfante operator 
Jz,bei (quark) used by Ji for the longitudinal component of the quark angular momentum, 
which has the nice property that it can be measured in Deeply-virtual Compton Scatter- 
ing reactions, does indeed represent the Z-component of the angular momentum carried 
by the quarks in a nucleon moving in the Z direction, and therefore, Ji's definition of the 
orbital angular momentum as the difference [J Zt bei (quark) — |AEjg^-], is fine as long as it is 
appreciated that this applies only to the components along the motion of the nucleon. 

E. The longitudinal component of the quark momentum and angular momentum 

We have argued that the canonical versions of the momentum and angular momentum 
operators should be regarded as the physically meaningful ones. Yet it is well known that xb, 
Bjorken-x, can be interpreted as the fraction of the Z component of the quark momentum 
inside a nucleon, in an infinite momentum frame where the nucleon is moving along the 
OZ axis, and that this corresponds, via the Operator Product Expansion, to the matrix 
element of the Bellinfante version of the momentum operators. At first sight this appears 
to contradict our assertion that it is the canonical version that should be regarded as the 
physically meaningful momentum operators. We shall here explain that there is, in fact, 
no contradiction in the special case of the longitudinal components of the momentum and 
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angular momentum. 

The gauge invariant expression for the unpolarized quark number density q(x) is usually 
written as 

= £ I ^e ixP+z ~ (P\^{-z-/2) 1 + W^{z-/2)\P)\ x>0 (113) 



2 J 2tt 

where | P ) corresponds to an unpolarized proton moving along the OZ axis i.e. 

P"=(P,0,0,P), (114) 

and where 

r z-/2 

W = W[-z~/2 , z~/2\ = V exp{i# / dz' A+(z'n) t a } (115) 

J-Z-/2 

is the Wilson line operator, a matrix in colour space, and where 

n= "^(1,0, 0,-1). (116) 
We are using the standard definition of the ± components of a vector i.e. 

v ± = -L(v + v z ). (117) 
The expression for the antiquark density q(x) is analogous to Eq. (11 13ft but with x < 0. 

After some manipulation one finds that 



x ^ x ) = 4^+ / ^e lxP+z ~ (P\{i>(z)[- < d + - ig A + (z)]} z= _ z - /2l + W^(z'/2) 

+ $(-z-/2) 7 + W{[S + - ig A + (z)] 4>{z)} z=z - /2 \ P)\ x>0 . (118) 



Integrating over x one has 



with 



<te [g(a;) + g(x) ] = P I ^(0) 1 + D + ^(0) | P ) (119) 



= ^+ - V+ - 2^A+(0). (120) 



Now from Eq. ( 1ST]) the quark part of t^(gG) is given by 

CwW = ^h^WW + ^ *)] - <T A (121) 
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where C q is the quark part of C q c given in Eq. ( 161]) . 
Then, since g ++ = we see that 



so that 



^ 1 dxx[g(x) + glx)] = ^ T ^(P|t+ + i (0) |P>- (123) 



Consider, now, the physical interpretation of the LHS of Eq. fl 1 2 3 [) in the parton model. The 
parton model is not synonymous with QCD. It is a picture, a manifestation, of QCD in the 
gauge A + = and it is in this gauge, and in an infinite momentum frame that x can be 
interpreted as the momentum fraction carried by a quark in the nucleon. But since A + = 
we have 

*L) + (gauged = 0) (124) 

so that for these particular components of the tensors there is no difference between the 
canonical and Bellinfante versions 

*J£»(0) = t+&(0) (gauged = 0). (125) 

Hence the fraction of longitudinal momentum carried by the quarks in an infinite momentum 
frame is given equally well by either the canonical or Belllinfante versions of the energy 
momentum tensor density. 

Let us turn now to the question of the angular momentum and, in particular, to Ji's 
relation of the quark angular momentum to the second moment of certaingeneralized parton 
distributions (GPDs) measurable in Deeply Vitual Compton Scattering [2j. In the standard 



notation (see e.g. the review of Diehl 25J) 



1 f dz~ 



2 J 2tt 
1 



(P'| V(-*"/2) ^Wip(z-/2) \P) 



2P+ 



iA 



2M 



£u{P , )a» p u{P) 



E(x,£,t) 



(126) 



where 



P = hp + P>) A = P'-P t = A 2 A + = -2£P^ 



(127) 
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and the spinors are normalized to uu = 2M. Putting P' — P i.e. A = t = £ = and 
comparing with Eq. fll 13j) one sees that 



if (x, 0,0) = q(x) (128) 

so that xH(x, 0, 0) can be interpreted as the density in x-space of the quark longitudinal 
momentum. 

Now consider the general expression for the off-diagonal nucleon matrix element of t^ u bel (0). 
The connection between these matrix elements and the angular momentum involves diver- 
gent integrals, which have to be treated carefully using wave packets, as was done correctly 



for arbitrary components of J for the first time by BLT 19j, and for this reason we shall use 



their notation for the scalar functions that appear in the matrix element of ^^(0). One 
has 



P', S' | I P > S ) = [u'YuP" + (/i u)]B qMl (A 2 )/2 



2M 



v) 



D 9;beZ (A 2 )/2-§ gifcei (A 2 )] 



u u 
2M 



J g , 6e/ (A 2 ) - H 9ibeZ (A 2 )](A^A^ - A 2 g^) + M 2 WL q , bel {A 2 )g^ 



(129) 



where 

u = u(P,S) u' = u(P',S'). (130) 

Note that the term M 2 M,g pu is only allowed because we are dealing with a non-conserved 
density. 

Repeating for the GPDs the analysis which led to Eq. ( 11231) and bearing in mind Eq. (I126p 
yields 

1 



2P+ 



u f/ y u] / dxxH(x, £, i) + 



iA P -/ +p 
— - u a H u 



1 



2M 

p',s'\KU q ) \P S )- 



dxxE(x, £, t) 



2{P+) 2 

From Eq. (11291) . remembering that g ++ = and that A + = — 2£P + , one obtains 

(P 1 , S ' I Kbei(°) \P^) = [u'j + u>P + ]p q ,U^) + e(G qMl (A 2 ) - H 9 , 6d (A 2 ))] 



(131) 



+ 



iA 



2M 



Wa +p uP^ 



[2§ q>bel (A 2 ) -© g , feei (A 2 ) -e{G qM i(A 2 ) -H,, beZ (A 2 ))]. (132) 
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Comparing with Eq. (113 ip . taking the limit A — > and writing V> q ^ei = ^q,bei 

(A 2 = 0) 

etc, one obtains 



J dxxH(x, 0,0) =3 qibd (133) 
J dxxE(x,0,0) = (2§ qMl -B qMl ) (134) 



and consequently 



J dxxH (x, 0, 0) + J dxxE(x, 0,0) =2§ q , bd . (135) 



Consider now the parton model interpretation of these expressions. Choosing the gauge 
A + = we have, as before, £^^ n (0) = ^(0), so that in Eqs. (ll33[ I134[ 1 1 35[) we may put 



"q, bel — >Wq, can = ^>q 



and \ bd = E> qiCan = S q . (136) 



For the case of a longitudinally polarized nucleoli moving at high speed in the Z direction 
BLT 



191 ] proved that S measures the Z-component of J. Hence Eq. (I135P can be written 
J dxx[H(x, 0, 0) + E(x, 0, 0)] = 2 J 2 (quark) (137) 
which is the relation first derived by Ji 

Note, however, that unlike the case of linear momentum, it is not obvious that x [H(x, 0, 0) + 
E(x, 0, 0)|| can be interpreted as the x-space density of J z (quark). Indeed, Burkardt and 
Hikmat [26] have shown, in a model, that J z (quark; x) calculated directly from the nucleon 
wave function disagrees with x [H(x, 0, 0) + E(x, 0, 0)]], whereas there is perfect agreement 
when integrated over x. 



F. Interacting particles: photons and gluons 

To a large extent the entire controversy concerning the assigning of angular momentum 
to quarks and gluons arose from the long established claim that one cannot split the angular 
momentum of a massless gauge particle into an orbital and spin part in a gauge-invariant 
way. The two standard expressions in the literature for the angular momentum for QED, 
the canonical and Bellinfante versions, are 

33 



Jcan = J d 3 X1p^J-f 5 1p + J (1 3 X^[X X (~lV)]V> 

+ J d 3 x(E x A) + J d 3 xE i (x x VA*) 

= S can {el) + L can {el) + 5 can (7) + £ C(m (7) (138) 



and 



'beZ 



+ J d 3 xx X (E X B) 

= S bel {el) + L bel {el) + J M ( 7 ) (139) 

In «7 can only the electron spin term is gauge invariant. In J be \ each of the three terms is 
gauge invariant, but the photon angular momentum is not split into a spin and orbital part. 
Insisting on being able to split the photon angular momentum into a spin and orbital part, 

n 

and on having each term gauge invariant, Chen et al [3] arrived at the following form 
J, 



d 3 x^jj 5 tp + J d 3 xip ] [xx (-iD puTe )}4) 
+ / d 3 x (E x A vhvs ) + / d 3 xE i (x x VA { 



chert 

' ^ ja - ,m - A phys ) + J a-xw\xx v/i physJ 
= S chen (el) + L chen (el) + Schen^) + L chen {i) (140) 

where D pure = V — ieA pure and the fields A pure and A phys were explained in Eqs. ([U13]) of 
the Introduction. 

Later, Wakamatsu [4] suggested a rearranged version of J c hen> which retains a gauge- 
invariant split between the spin and orbital angular momentum of the photon 



Jwak = J d 3 xip^'j 5 ip + J d 3 xi) ] [x x {-iD)\ij} 

+ J d 3 x (E x A phys ) + [ J d 3 xE\x x VA phys ) + J d 3 x^(x x eA phys )ip) 

= S wak (el) + L wak (el) + S wak (j) + L wak (j) (141) 

In this version the very last term J d 3 xif} J *(x x eA p h ys )ip has been shifted from Chen et al's 
electron orbital term to the photon's orbital angular momentum. We have already com- 
mented that one could do such a rearrangement in an infinite number of ways by shifting 
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some arbitrary fraction of this term. 

All of the above comments hold equally well for the case of QCD. 

As we have stressed, there is absolutely no need to have gauge-invariant operators so long as 
their physical matrix elements are gauge invariant, as is the case for the canonical version 
of the electron spin, the electron orbital angular momentum and the photon's total angular 
momentum. Moreover we have insisted that an angular momentum operator should generate 
rotations, at least in the restricted "minimal" sense defined in Eq. (1931) . Only the canonical 
choice satisfies this requirement. We conclude, in agreement with the paper of Jaffe and 
Manohar [7], that it is the terms in the canonical form J can which should be interpreted as 
corresponding to the angular momentum of the electron and photon respectively. Of course 
this leaves open the issue of splitting the photon angular momentum into spin and orbital 
parts. This we shall discuss in the next section. 



VII. THE SPIN OF THE PHOTON AND THE GLUON 

As has been emphasized for more than half a century it is true that the canonical photon 
or gluon spin terms, as a whole, are not gauge invariant. This we regard as an inevitable 
feature of a gauge theory and it has not been the cause of any problems in the description 
and calculation of physical processes involving photons, and more recently, gluons. However, 
the projection of the spin terms onto the direction of the photon's or gluon's momentum 
i.e. the photon and gluon helicity, is gauge invariant and it is this quantity which can 
be measured and, as we shall show, is measured in deep inelastic scattering on atoms or 
nucleons respectively. 

A. QED 

Consider the expression for S ca n(l) in Eq. (11381) . which can be written as 

SL(l) = \^S^. (142) 

In accordance with the analogue of Eq. (11031) we may study its expectation value between 
states of definite mass by replacing the fields by their "in-field" versions Thus we may utilize 
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Sg = jd 3 x [FL°(x)4 n (x) - Fi°(x)AUx)}. (143) 
We shall show that the matrix element of the "in-field" helicity operator 

^n=[S k can {l)P k l\P\] m (144) 

taken between arbitrary physical states of a photon is gauge invariant and then relate its 
matrix elements to the QED analogue of the polarized gluon density Ag(x). 
Consider the action of S\l on a physical photon state of momentum k and polarization 
vector e(k , /) corresponding to polarization along a transverse direction I : 

\k,l) = a\k,l) | vac). (145) 
Provided the operators are normal ordered we have, 

^|fc,/) = [^ ; at(fc,/)]|vac). (146) 

Then if i and j correspond to directions perpendicular to k, expressing a^(k , I) in terms 
of the fields as in Section 14.4 of [9|, and using the equal time commutators (permitted 
because we are dealing with the "in- field" momentum and angular momentum), gives 

[SU , at(fe , I)] = i{e\k , I) a\k,j) - e*(fc , I) o+(fc , i)}. (147) 

Let us first check that acting on a helicity state, %i n , as given by Eqs. (11441 11421) . yields 
the correct result when using Eqs. (I146[ I147p . For simplicity take the OZ axis along k. 
Then 

U m = [S cantZ { 1 )] m = S}Z (148) 
and the helicity states are, for A = ±1, 

\kz, A) = ^{\kz, l) + iX\kz, 2)}. (149) 
v2 

Using the fact that e^(kz , 1) = 8j% , eP(kz , 2) = 5 j2 one finds that indeed 

n in \kz, X) = X\kz, X). (150) 
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To show the gauge invariance of we consider its action on a general physical photon 
state 

|$) = ! d 3 kJ2Mk)\k,l) (151) 

l±k 

where the sum over / refers to directions perpendicular to k. Then 

U in \®)= 1 - j d 3 kJ2Mk)^e n3 SH\kJ). (152) 



l±k 

Since i, j and I refer to directions orthogonal to k we may use the results Eqs. (1146} I147j) to 
obtain 

n in | $ ) = l - / d 3 k Hk) ttt eni {e\k , I) a\k,j) - e j (k ,l)a\k,i)}\ vac ). (153) 

J l±k ' ' 

The creation operators in Eq. ( 11531) refer to polarization directions orthogonal to k, and are 
thus unaffected by gauge transformations. Hence the most general matrix element of Hj n is 
gauge invariant. 



For QED we now introduce A7(x), t he g auge-invariant ana 



ogue of Ag(x), based on the 
expression for Ag(x) given by Manohar 27j and used by Jaffe ^1 i.e. 

A 7 (x) = j dCe-^- p+ (P, Sl |^ a (D/(r,0)F+ im (0) |P, S L ) + (x -> -x) 

(154) 

where I(£~~,0) is the Wilson line integral, and |P, Sl) is a longitudinally polarized, fast 
moving state. The axes are chosen so that P = (0, 0, P) and we have used Eq. (I103j) to 
replace the fields by their "in-field" versions. 

Since the expression in Eq. (I154p is gau ge invariant we may evaluate it in the gauge 
A + = 0. Then following the argument in 28| and integrating over x, we obtain 



A7 = / dxAj(x) 



' (P, S L \Fl n + (0)Al(0) - F^(0)A] n (0) |P, S L ). (155) 



2P^ 

Consideration of the possible tensorial structure for the matrix elements indicates that 
in leading twist 



7 Note that there is a typographical error in the expression for Ag(x) in these papers: G a + (0) should be 
G+JO) 
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(P, S L \F l0 (0)\P, S L ) = (P, S L \F i3 (0)\P, S L ) (156) 
so that in leading twist 

A7 = ^ (P, S L \F? n °(0)Al(0) - ^°(0)A 1 (0) m \P, S L ). (157) 
Now using Eq. (11431) one sees that 

(P f , S L | S% |P, S L ) = (2vr) 3 5(P' - P) (P', S L \Fl n °(0)Al(0) - F™(0)A\0) in |P, S L ). 

(158) 

Hence 

A7_ 2E(2tt)3 tf(P'-P) • (i59j 
But the denominator is just the norm of the state |P, Sl) so that A7 indeed measures 
the expectation value of the photon helicity operator. 

B. QCD 

Because we may use the "in-fields" to study the matrix elements of the gluon helicity 
between arbitrary states of a nucleon, there is no essential difference from the photon case. 
The expression Eq. (11431) for S 1 ^ is simply altered by adding a colour label to the fields and 
summing over it. Similarly the expression for Ag(x) and Ag are obtained from Eqs. (I154p 
and (I157p by adding colour labels and summing over them. 

Thus Ag indeed measures the expectation value of the gluon helicity in a nucleon. 

VIII. CONCLUSIONS 

We have argued that there is no need to insist that the operators appearing in expres- 
sions for the momentum and angular momentum of the constituents of an interacting system 
should be gauge invariant, provided that the physical matrix elements of these operators 
are gauge invariant. We have also suggested that the expressions given by Chen et al and 
Wakamatsu for the momentum and angular momentum operators of quarks and gluons are 
somewhat arbitrary and do not satisfy the fundamental requirement that these operators 
should generate the relevant infinitesmal symmetry transformations specified in Eqs. fl92| 
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[931) . Demanding that the conditions Eqs. (|92| 193]) be satisfied leads to the conclusion that 
the canonical expressions for the momentum and angular momentum operators are the 
correct and physically meaningful ones. 

It is then an inescapable fact that the photon and gluon angular momentum operators 
cannot, in general, be split in a gauge- invariant way into a spin and orbital part. However, 
as discussed in detail, the projection of the photon and gluon spin onto their direction of 
motion i.e. their helicity, is gauge-invariant and is measured in deep inelastic scattering on 
atoms or nucleons respectively. 

Although Ji's expressions for the quark and gluon angular momenta, which are the Bellinfante 
versions, do not conform to the above conditions and thus should not be considered as 
measuring arbitrary components of the quark and gluon momenta and angular momenta, 
nonetheless, it turns out that the expectation value of the Bellinfante operator J Zy bei (quark) 
used by Ji for the longitudinal component of the quark angular momentum, which has 
the nice property that it can be measured in Deeply-virtual Compton Scattering reactions, 
does indeed represent the Z-component of the angular momentum carried by the quarks 
in a nucleon moving in the Z direction, and therefore, Ji's definition of the orbital angular 
momentum as the difference [J z> bei (quark) — |ASjg^-], is fine as long as it is appreciated 
that this applies only to the components along the motion of the nucleon. 
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